Single photon nonlinearities using arrays of cold polar molecules 
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We model single photon nonlinearities resulting from the dipole-dipole interactions of cold polar 
molecules. We propose utilizing "dark state polaritons" to effectively couple photon and molecular 
states; through this framework, coherent control of the nonlinearity can be expressed and potentially 
used in an optical quantum computation architecture. Due to the dipole-dipole interaction the 
photons pick up a measurable nonlinear phase even in the single photon regime. A manifold of 
protected symmetric eigenstates is used as basis. Depending on the implementation, major sources 
of decoherence result from non-symmetric interactions and phonon dispersion. We discuss the 
strength of the nonlinearity per photon and the feasibility of this system. 



I. INTRODUCTION 

Coherent control of optical nonlinearities at the sin- 
gle photon level is a burgeoning topic in quantum optics 
research. Utilizing state-preserving techniques, it is sug- 
gested that one can implement two-qubit quantum logic 
gates in a feasibly robust optical quantum computational 
framework [U |2] . 

Cold polar molecules are excellent candidates as a me- 
diating medium due to their field-dependent intermolec- 
ular interaction properties [3l HI |5l |6j |7]. They have 
been suggested for quantum computation architectures 
since they embody advantages of both neutral atoms and 
trapped ions, viz. long coherence times and strong inter- 
actions, respectively. 

Advances in preparation (cooling and trapping) of 
molecular ensembles in their electronic, vibrational, and 
rotational ground states [8] would allow for single state 
manipulation in a characteristically rich level structure. 
Notably, recent work by Biichler et al. |3] predicts novel, 
controllable superfluid and crystalline phase transitions 
from dipolar gases. The latter could suppress dephasing 
from short range collisions in high density traps. The 
anisotropic and long-range form of the dipole-dipole in- 
teraction is responsible for the bulk of advances in con- 
trolling molecular samples |5j. 

In this paper, we investigate cold polar molecular gases 
in one- and two-dimensional arrays. We describe single- 
photon nonlinearities resulting from the intermolecular 
dipole-dipole interaction. We apply "slow" and "stored" 
light methodology for coherent state transfer. Next we 
calculate the resultant nonlinear phase in the context of 
collective excitations in an optically thick media. 

Here we are primarily concerned with exploring fea- 
sibility of coherent control over the resultant nonlinear 
phase evolution of intermolecular dipole-dipole interac- 
tions. Familiar implementations for the system under dis- 
cussion include stripline cavities, optical lattices, Wigner 
crystals, hollow fibers, or molecules on surfaces. Then, 
we investigate the most significant decoherence effects for 



implementation in a trap architecture or in a crystalline 
phase. 

II. SINGLE-PHOTON NONLINEARITY 

Photons do not interact. However, effective interaction 
can be achieved by utilizing state-preserving light-matter 
couplings to nonlinear media, wherein matter-matter in- 
teractions effectuate photon-photon interactions without 
destroying the state information of the incident coherent 
fields. 

The proposed mechanism is as follows: Photons are 
efficiently and coherently coupled into the molecular 
medium in the form of "slow-light polaritons." The 
molecule part of these light-molecule coupled excitations 
is "switched" from the zero-dipole rotational ground 
state into a high-dipole rotational superposition state. 
The resulting dipole-dipole interaction therefore adds a 
nonlinear phase to the polaritons that the photons re- 
tain on exiting the medium. We note that the nonlinear 
phase is thus proportional to the interaction time inside 
the medium and thus to the propagation time of the po- 
laritons. Therefore control over the phase is exercised by 
manipulating the propagation velocity of the light in the 
medium. 

Electromagnetically induced transparency (EIT) - 
based slow light polaritons [IQl HI] are collective states of 
matter-light superposition that can be achieved by using 
a A-type system. Polaritons are the coupled exchanges of 
the signal field Q.s and the superposed \g) and |e) ground 
states (see Fig.jlJ. \g) and |e) would typically be the 
\J,Mj) rotational states of the ground state molecules, 
where Mj is the projection of J on the ^-axis. The cou- 
pling field Vic controls the slow group velocity of the po- 
laritons. 

The interacting states \g) and |e) in our system are 
neighboring rotational levels of dipolar molecules. Since 
dipolar interactions exchange virtual photons, interact- 
ing states must have opposite parity. This can be accom- 
plished via an expanded A-type system with a strong 
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FIG. 1: Level scheme for utilizing slow-light polaritoris. |e) 
and I a) can be coupled via a two-photon transition, or |e) can 
be a mixed parity state. 



Raman transition involving a two-photon transition as 
f2c, or alternatively the use of mixed-parity states for the 
excited state \a). For simplicity we refer to the whole 
molecular system as an effective two-level system, con- 
sisting of |e) and \g), as is usually done in the context of 
slow- light polaritons \W . We adopt a natural shorthand, 
\gi) and \ei) respectively for the ground and excited state 
of the ith molecule. 

It is convenient to introduce collective states denoted 
as \j,m). These states are eigenstates of the collective 



operators and Jz , where — }^ ^ . af , < 
and trf are Pauli operators acting on the i*^ molecule: 
= \9^) (e»| + |ei) {gr\, o\ = i{\g^) {ei\ - \ei) {gi\), &l = 
|ei) {e.i\ — \gi) {gi\. | j, to) states satisfy the eigenvalue rela- 
tions to) = j(j + 1) |j, to) and Jz |j, m) = m \j, to), 
with j — N/2, . . . , and —j < m < j. 

Among these states we are particular interested in the 
fully symmetric Dicke-like states which lie on the surface 
of the Bloch sphere with maximal radius j — N/2 and are 
totally symmetric, i.e. invariant with respect to particle 
permutations. We denote them as \n) = \N/2, —N/2 + n) 
to emphasize that they correspond to n-photon collective 
excitation. The corresponding n = 0,1,2 are explicitly 
given by 
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In our case, cold polar molecular gases e.g., SrO [5] or 
CaF [12] constitute the nonlinear medium. The non- 
linearity is expressed through dipole-dipole interactions. 
Imagine first the ideal case when such interactions be- 
tween molecules generate an effective Hamiltonian of the 



Vdd = xJz 



(2) 



Since the states \n) are eigenstates of Vdd, the collective 
dynamics can be fully accounted for by their phase evo- 
lution 9„'. 



he,,{t) = (n| Vddt \n) = x(^/2 - nft 



(3) 



In general to characterize the medium nonlinearity we 
would have to include photon coupling states with n > 
2, but in the scheme of optical quantum computation 
it is sufficient to implement two-qubit controlled phase 
operation - provided appropriate single qubit gates. The 
latter can be realized if the accumulated non-linear phase 
B(t) acquired by the polaritons equals tt at the time when 
they exit the medium. Q{t) is the difference between the 
phase picked up by two concurrent excitations and the 
sum of the phases that each individual excitation would 
independently pick up in the absence of the other. It is 
defined as 

e(t) = (02(i)-0o(t))-2(0i(t)-eo(t)) - 92-261+90. (4) 

In other words, 0(t) quantifies the departure from a lin- 
ear regime, that is, one in which the interaction between 
the two excitations is absent and therefore their indi- 
vidual phases simply sum up. From Eq. ^ the latter 
condition is satisfied if 2xt.„ = fnr. Therefore in this 
ideal case establishing a deterministic controlled phase 
operation only requires coherent control of the propaga- 
tion and/or storage time of the polariton in the dipolar 
medium. This corresponds to manipulating the control 
fields that establish the conditions for the "slow light" 
propagation. 

However, dipolar interactions do not generate a Hamil- 
tonian of the type described by Eq. Q and instead the 
dipole-dipole interaction is given by 
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with fii the dipole moment of the molecule at site r°. 
Here we have assumed that the molecules are at fixed 
positions determined for example by a superimposed ex- 
ternal optical lattice potential. 

For the simplest situation when both \g) and |e) states 
are pure rotational states and have zero dipole moment 
(j,gg = /.lee — 0; the dipole-dipole interaction is governed 
by the \g) <-> |e) transition dipole moment which can be 
formally written as: 
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with /inf, = fiQ {a\er\b). Assuming the interacting 
dipoles are aligned in parallel, which is possible in ID and 
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2D geometries, this leads to {p,i ■ r°) = 0, and neglect- 
ing counter-rotating terms d'fa^, <y^(J^ the interaction 
potential becomes 



dd ■ 



r.0|3 



(7) 



Vdd is not SU(2) symmetric and consequently the col- 
lective Dicke states are not eigenstates of it. Exceptions 
are |0) and |1) states which do remain eigestates of Vdd- 
This implies that the dynamical evolution of |n) for n > 2 
not only acquires a time-dependent phase, but in addi- 
tion transitions to other states outside j = N/2 will take 
place. These transitions will affect the implementation 
of the phase gate which relies on remaining on the Dicke 
manifold. 

Ignoring for the moment the "leakage" outside the 
Dicke states and focussing only the projection V of Vdd 
on the Dicke manifold, which is given by [13 
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and a the lattice constant of the molecular array. Now 
one can estimate the propagation time required for im- 
plementing the phase gate as t^^ = hTT/(2xcs)- The re- 
sulting expression gives Xeff oc Yjv- It clearly shows that 
there is an optimization to undertake regarding the num- 
ber of molecules in the array: On one hand, there must 
be enough molecules to create sufficient optical depth to 
couple-in the polaritons |14l . On the other hand, in order 



to maximize the nonlinearity in Eq. ( 10 1 , less molecules 
are better. 

In a one dimensional (ID) molecular array, one can 
analytically evaluate the nonlinear phase factor O from 
Eq. It is given by : 



h{N-i) 



(10) 



where C[3] = limAr^oo Y.l=i ^ « ^ « 1.2. 

Given aforementioned assumptions about dipole align- 
ment, for the two-dimensional square lattice [T7] only a 
change of lattice vectors — yiBy + Zie^ is required. In 
this case and assuming periodic boundary conditions one 
obtains 

9(215)^29(1^) (11) 



III. DECOHERENCE 

a. Decay Out of Symmetric manifolds 

Decay out of symmetric manifolds and phonon-like ef- 
fects in the Wigner crystal implementations can be sig- 




FIG. 2: (Color Online): (a) Nonlinear phase as a function of 
time for a ID array with = 36 molecules (blue solid line) 
and iV = 81 (red dashed line). Here U(N = 36) = fi7r/(2xoff) 
is the expected phase gate time for the A'' = 36 system cal- 
culated by projecting the dipole Hamiltonian (Eq. (|5|) onto 
the Dicke manifold. It is obvious from this graph, that the 
projected Q — yr-phase time {t = f-^, solid black line for 
N = 36, broken for A*' = 81) deviates strongly from the exact 
one {t « 3.5 t,r). The main reason for this deviation is the 
importance of transition processes out of the Dicke manifold 
as confirmed in panel (b) where we plot the fidelity to stay in 
the state |2). 



nificant. We will discuss symmetric manifolds at present, 
leaving the phonon-like decoherence effects to a later sec- 
tion. 

As mentioned in previous session, Dicke states are a 
good basis only if the relevant Hamiltonian is spherically 
symmetric (SU(2) symmetric). This is not the case for 
Vdd, and in particular the state |2) will decay during the 
time evolution inducing decoherence. 

To estimate the decay probability from an initial Dicke 
eigenstate during the dynamical evolution we calculate 
the fideHty F{t) 



m = 



m) |2) 

(^(0) |2) 



(12) 



where |V'(^)) = ^f^dd |V'(i))i is the time evolving state 
under Vdd and |^(0)) = ^2(0) |2) + CM |1) + Co(0) |0). 
This quantity is plotted for two different ID sample sizes 
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m Fig.[2| In the same figure we also show the non-linear 
phase Q accumulated by the evolving state. We numeri- 
cally evaluated it as 



cos[e(t)/2] 



(13) 



where C„(i) = ((V'(^) l^))/(V'(0) 1*^) are the projections 
of the evolving state into the corresponding Dicke states. 
The relevance of decoherence effects and the departure of 
the pure phase accumulation can be clearly observed in 
Fig.|2] The plot shows not only a distorted evolution of 
the non-linear phase but also a different time dynamics 
since the non-linear phase approaches tt at a time very 
different from the expected (See Eq. (10)) . 

One could effectively remove the mixture of j manifolds 
and improve the fidelity of the phase gate by the addition 
of an external electric field. This procedure, which we will 
describe in the following section, generally establishes a 
Many-body Protected Manifold (MPM) [fsl which helps 
to eliminate or mitigate decoherence effects. 



b. Phonon-like effects 

In the presence of an external DC field, which induce 
repulsive dipole-dipole interactions in the ground and ex- 
cited states, (^gg, ^ee 7^ 0), molecules can assemble them- 
selves in a Wigner crystal. Attractive interactions along 
the remaining directions can be suppressed by a strong 
transverse confinement [9]. In a Wigner crystal imple- 
mentation, there exists another notable decoherence ef- 
fect, which can be analyzed using a phonon formalism, 
cf. [T5]. In a realistic crystalline phase, the molecules are 
not fixed frozen. Phonon-like effects will add to the decay 
described in the previous section as their energy provides 
a coupling between the symmetric and non-symmetric 
states. Details of decoherence due to phonons will also 
be treated in subsequent sections. 



c. Finite pulse effects 

In the previous analysis we have assumed that Dicke 
states are the result of the coherent light-molecule inter- 
actions. However finite pulse effects can introduce in- 
homogeneity and can lead to non-zero initial population 
of states out-side the Dicke manifold. The latter will 
cause additional decoherence and will degrade the phase 
gate. The corrections can be quantitatively understood 
by noticing that while slow light polaritons have linear 
dispersion in a linear medium, the nonlinear interaction 
adds a dispersion relation: 



Eik) = hu;^ = «E 1^0 013 • ^°)- 



For ID, Eq. (14 1 can be rewritten as 



ID 



K ((-3 -I- 2 In ka) {kaf + 0{kf 
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The non-linear terms in k present in the ID dispersion 
relation will degrade the phase gate. They however can 
be mitigated by using long pulses or a (ring) cavity, where 
fc = 0. 

In 2D on the contrary the low energy excitations scale 

as 



hjjj 



2D 



3.27K|ka| 



(16) 



showing that in the 2D case, at least in the long-wave 
limit, the spectrum remains linear and thus decoherence 
due to finite pulse effects becomes less important. 



IV. MANYBODY PROTECTED MANIFOLD 
(MPM) 

As the next step, we include a tunable DC electric field 
which thereby induces a dipole moment in the ground and 
excited states of the molecule {ngg,^ee 0)- This effect 
then augments the dipole-dipole interaction among our 
collective Dicke-like states in a way that enables pertur- 
bative treatment of the non-spherically symmetric part 
of the interaction, thus reinstating |j, to) as good eigen- 
states for the system: 



V 



Hh + Hi 



(17) 



where ^ 



iMegl 2 (/^'^'^ f^gg) 



Here Hh is the spherically symmetric (Heisenberg) part 
of the Hamiltonian V, and Hj is the non-symmetric 
(Ising) part. 

If at t = an initial state is prepared within the 
j = N/2 manifold, a perturbative analysis predicts that 
for times t such that Kt/H < £,/k, Hh confines the dy- 
namics to the Dicke manifold and transitions outside it 
can be neglected. In other words the Dicke manifold 
becomes protected by the many-body interactions and 
only the projection of Hi on it, which corresponds to 
VHj = —XcsJz + const, with 

Xeff = f Xoff becomes effective. As a consequence Hj 
acts as the desired ideal "phase gate" Hamiltonian. 

The relative strength of the Hh and Hj parts of V can 
be manipulated to find values of ^ /k such that MPM pro- 
tection is maximized. For example, SrO has a ground 
state with a magnetic moment of 8.89D. Upon selecting 
ground and excited rotational states with opposing par- 
ity, the appropriate values of ^ and k are then obtained 
by diagonalizing the Stark Hamiltonian for variable elec- 
tric fields [13. In Fig. [3] we show the ratio of ^/k for 
varying DC field strength between two rotational levels 
of SrO. For induced dipole transitions, the biasing elec- 
tric DC field E depends on both the rotational constant 
B and the ground state dipole moment /xq. 



5 



y 



-OS 

-1 




0.5 1 1-5 

(a) 



FIG. 3: C/k for varying E[B/fio] in SrO. (a) \gi) = 
|JsrO,Mj,s.o) = |0,0), and |e,) = |1,0),. (b) \g,) = |1,0), 
and |ei) = |2, 0)^. B is the rotational constant of the molecule, 
and fj,o is the maximum ground state dipole-moment, in our 
case assumed to be ~ iMegl- 



The application of an electric field changes the original 
bare states to dressed states \g) and |e) which are linear 
superpositions of the bare states. It has to be noted that 
the addition of a DC field, leading to MPM protection, 
causes the nonlinearity to be reduced by a factor of ^ 



Decay Out of Manybody Protected manifold 
into other manifolds 



In this section we study what occurs if at time t ^ 
we prepare the system in the j — N/2 subspace and let 
the system evolve in time in the presence of MPM. This 
can be written as 



= e-^^"(*)*e-«"(*) r/2,~72 + ") (18) 



Here we estimate the decay magnitude 7„. Our objec- 
tives are to keep 0„ large while minimizing 7„. As the 
states with and 1 excitations are eigenstates of the sym- 
metric manifold there is no decay out of the manifolds 
and therefore 70 = and 71 = 0. 

Using first-order perturbation theory we can write 



-72 (t)t 



-si: 

fc,fc'>0 



dTMk,k'e 



(19) 



The quantities Mk,k' = (^/2, -^/2 + 2| i// |V'fc,fc'>o) 
are the transition matrix elements to states with j — 
N/2 — 2 which are the only ones which couple to 
1 2) according to the Wigner-Eckart theorem. To a 
good approximation they are given by: |V'fe,fc') = 

, ^ : V- / ■ e'^''"''°J''"'°'"^^V~g~|0) and their corre- 
^N(N-l) J * ' ' 

sponding excitation energies by h{ujk + lju'), with LOk 

given by Eq. (14). hk, hk' are discrete quasi- momenta. 



which for a 2D square lattice with lattice spacing a can 



be written as k 



/N 



i,j = 0,...N-1. Note that 



the sum over k, k' in Aik,k' excludes the state k = k' = 
since ■00,0 is just |2). 

After some algebra, one can show that A4k,k' — 

-^Fk5k,—k' where Fk is the Fourier series of — r^j^'^. 



Fk = 



r^l cos(fc • r"). Replacing the 



latter equation in Eq. ( 19 1 yields the following expression 
for the decay rate 
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To get a general idea on the decay rate behavior we 
first use the Fermi- Golden Rule to estimate the decay 
rate in the thermodynamic limit, N 00 and then we 
compare this predictions with numerical studies for finite 
size systems. According to the Fermi-Golden rule, at 
long times the decay probability evolves linear with time 
as 72 (t) = 
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The latter relation yields that r2tTr diverges in ID as 



(21) 
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implying the break down of the Fermi-Golden rule ap- 
proximation and emphasizing the issue that in ID non- 
symmetric decoherence effects are crucial in the large N 
limit. 

In 2D, the situation is better due to the linear de- 
pendence of the long wave excitations with k and the 
extra-factor of k in the density of states. This yields that 



r2i^ CX 



d{ka)kaS(ka) — > 



and 



F'^iU) - 1, 



(22) 



Consequently as long as ^ < k (which is required for 
the validity of our perturbative treatment) and neglect- 
ing other decoherence effects during the time evolution 
(which grows linearly with N) the Fermi-Golden rule pre- 
dicts a robust phase gate in 2D. 

To validate this predictions we solve the exact many- 
body dynamics numerically by evolving a system initially 
prepared in the Dicke state with n ~ 2 under the effective 
Hamiltonian V, and compute the fidelity, J-' of remaining 
in a Dicke state: 



m= r/2,-^/2+2|vI/(t)) 



= p-72(t)* 



(23) 



In Fig.|4] we show the ID dynamics using the param- 
eters £,/k = 0.05 and iV = 36 and 81 and plot both Q 
and F{t). Note in the presence of the MPM, the time 
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when the phase gate is implemented is close to the ex- 
pected time t^r- Fig. [4] confirms the prediction that in ID 
the fidelity decreases with increasing TV. We find that 
for moderate N the ID decay probability increases as 

e 



which is obtained via a fit. This relation 



implies that in order to implement a robust gate 



100 



j2 ^ ^1.62 



(24) 



However, by choosing ^/k small we pay the price of hav- 
ing slower dynamics and therefore we make the system 
more vulnerable to other type of losses. 

Fig.[5]emphasizes the gain in fidelity obtained by going 
from ID to 2D. With the same number of molecules and 
even a much larger ^/k = 1 the fidelity is much better 
than in ID. By numerically evaluating the maximum de- 
cay probability we find it behaves as (x ^N^'^-^^. The 
decrease in F with increasing N is in agreement with the 
Fermi-Golden Rule approximation. 



V. WIGNER CRYSTAL AND PHONONS 

In dense low temperature systems with sufficiently 
strong fixed DC electric fields, [3], one can realize a self 
ensemble molecular crystal or Wigner crystal. In this 
crystalline phase |^gg| > and the molecules are local- 
ized at their classical equilibrium positions, r°. The lat- 
ter form a linear chain in ID or a triangular lattice in 
2D, with lattice spacing a. The formation of a Wigner 
crystal is fundamentally determined by the dimensionless 
parameter 



potential energy 
kinetic energy 



dd 



(25) 



for molecules of mass m for a given density p — l/{a)^ . 
For P ^ 1 the dipolar repulsion wins over kinetic en- 
ergy stabilizing the crystalline phase. In contrast to the 
case where the the localization of the molecules is en- 
forced by an external potential such as an optical lat- 
tice, in the self-assembled crystalline phase, molecules 
are not frozen and they can exhibit collective oscillations 
(phonons) about their equilibrium positions. These os- 
cillations can be described by rewriting the position op- 
erators as ri = r'^ + Xi and expanding V in powers of 
Xi. This procedure yields three terms: the fixed position 
dipolar Hamiltonian (Eq. ([T8|) described in the prior sec- 
tions, the phonon Hamiltonian which is gapless (see Ap- 
pendix A) and a phonon-polariton interaction Hamilto- 
nian Vphon-po given by 



Vphon-po — — 2k ^ ^ GijCJi ■ aj — 2^ ^ ^ GijCj^a^ 



8B, 



(26) 
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FIG. 4: (Color Online): a) Nonlinear phase as a function of 
time for a ID array with = 36 molecules (blue solid line) 
and A'^ = 81 (red dashed line) in the presence of an external 
DC field, ^/tt — 0.05. The latter is used to implement the 
MPM. Here t^(7V = 36) = ft7r/(2xeff) is the expected phase 
gate time from our perturbative analysis for the A*' = 36 sys- 
tem (indicated by a solid grid line). The corresponding time 
for the A'' = 81 system is indicated by the dashed vertical 
line. For the two cases the actual time at which the phase 
gate is accomplished, |cos(0/2)| = is close to the calcu- 
lated indicating the validity of the perturbative analysis 
specially for A'^ = 36. The deviations can be accounted for by 
higher order corrections in perturbation theory. The fidelity 
of remaining in the |2) state is shown in panel b. The latter 
decreases as either the ratio f/« or A' increases, consistently 
with the Fermi-Golden rule predictions. 



where 



-3a- 



(xi-Xj).(rP-r9) 



4|r? 



and Bo = {pl^ - /igg)/(87ra^eo). Vphon-po is not spheri- 
cally symmetric, Dicke states are not eigenstates of it and 
consequently Vphon-po will induce transition outside the 
Dicke manifold even for the one-photon excitation state, 
|1). This transitions can degrade the phase gate signif- 
icantly since they are not suppressed by the MPM due 
to the gapless nature of the phonon spectrum. If at time 
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t/t7r(N=36) 
(b) 

FIG. 5: (Color Online): a) Nonlinear phase as a function of 
time for a 2D array with N = 36 molecules (blue solid line) 
and TV = 81 (red dashed line). The solid and dashed vertical 
lines indicated the expected phase gate time according to our 
perturbative analysis. The ratio of ^/k = 1 is outside the 
perturbative regime however the dynamics exhibits almost the 
behavior expected from an ideal evolution. The fidelity to 
remain in the |2) state is shown in panel b. Note that in 2D 
the fidelity improves as TV is increased. 



t = we prepare our state in the J = TV/2 manifold, the 
projection of the evolving state on the n— Dicke state 
can be written as 



Here N{ujxOi)) = l/(e'*"^('')/(feBT) _ jg ^^le thermal 
occupation number for phonons with the phonon spec- 
trum ujxik) (see Fig.|8]and Fig.|9]). See Eq. ^ for the 
definition of g\{k). For two photon excitations it can be 
shown in a similar way that 



l2ph ~ 27 



Iph- 



(29) 



Detailed derivations are included in App. B. 

In Ref. [15] analytical expressions for the decay rates 
in the thermodynamic limit were derived using the Fermi 
Golden Rule. In this limit the decoherence induced by 
phonons was shown only to be relevant in ID and pro- 
portional to the temperature. For our finite number of 
molecules Fermi Golden rule results are only crude ap- 
proximations and for a quantitative treatment we instead 
perform numerical calculations which are summarized in 
Fig. |6] and Fig. [7] There we plot the time evolution of the 
decay probability for two systems with different TV. We 
find a general tendency of the maximum decay probabil- 
ity to grow with increasing TV specially for ID. Also for 
fixed TV, (3 and k we observe a linear dependence on the 
temperature in agreement with the Fermi Golden Rule 
predictions. 
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t/t7r(N=36) 



3.0 



FIG. 6: (Color online) l-Fi(t), where l-Fi{t) / {i,+'iBof / ^ 
is the probability of decay of a single excitation for a ID 
system with TV = 36(blue), TV = 81 (red-dashed) and ^, Bo and 
13 are system dependant parameters. The solid and dashed 
vertical lines are at the time where the perturbative treatment 
predicts the implementation of the phase gate for the TV = 36 
and 81 systems respectively. 



|vl,„(i)) = e-^*e-«"W + n) . (27) 

and the decoherence rates ^ipu and ^2ph approximately 
calculated using perturbation theory |15j . This proce- 
dure yields an expression for ^ipu given by: 

+N{ujx{k))5{ujx{k) + Uk)]. 

(28) 



VI. CONCLUSION 

In this paper, we explore the feasibility of utilizing cold 
polar molecules in ID and 2D optical lattices for coher- 
ently controlled nonlinear optics. We report a controlled 
TT phase gate time that increases proportionally to the 
number of interacting molecules, but also note better fi- 
delity in 2D systems for reasonable system parameters 
and external field strengths. We address the role of non- 
symmetric interactions, one of the major decoherence ef- 
fects, and demonstrate the enhancement of phase gate 
fidelity when an MPM is created by applying external 
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t/t;r(N=8l2) 

0. 0.2 0.4 0.6 0.8 1. 
0.25iH 1 1 1 1 h 




t/t;r(N=202) 



FIG. 7: (Color Online) 1 - Fi{t), where 1 - Fi(t)/{^ + 
4_Bo)^/\//? is the probability of decay of a single excitation 
for a 2D system with = 20^(blue), N = 81^ (red-dashed) 
and ^, -Bo and /3 are system dependant parameters. In this 
case we had to go to larger system sizes to show the decrease 
in fidelity with increasing A^. The top (bottom) axis scales 
are for N = 81^(20^). 

electric fields. For self assembled crystalline samples, we 
also have explored phonon-induced decoherence. Since 
we find that at low temperature the most relevant deco- 



herence effects in the Wigner crystal arrays are caused 
by long wave phonon excitations, spin echo techniques 
could help to reduce them. 
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Appendix A 

In a self assembled ensemble molecular crystal, 
molecules are not longer completely frozen at the classi- 
cal equilibrium positions, , which form a linear chain 
in ID or a triangular lattice in 2D, with lattice spacing a. 
Instead they exhibit collective oscillations (phonons) in 
the crystal. If we expand the total Hamiltonian, Eq. ([T]) 
around the equilibrium positions: rj = rP -I- Xi and keep 
terms up to quadratic order Xi one obtains the following 
expression 



H — Hphon + V+ Vphon-po (30) 

rr _ yPf , ^ 5[(x.-xj).n°/-[(x.-xj)]^ 

tiphon - 2^ 2to + 167reoa3 2^ \v9 - r9|5 ^■^'> 

= '^f«^xiQ)al^x^q.\ (32) 



where ciq^x is the annihilation operator for phonons can write 
with quasi-momentum q and frequency u!\{q) and is 
a unit vector along r? — r9 . In 2D the index A labels the 
two different phonon branches. In general 



t ft' 



ux{q)=^h{q) (33) e-^^Mt)^ i _ ^ j dt' j dr l^'^.iM'^ 



Jo 



q,X,k>0 

The phonon modes in the dipolar crystal are acoustic [{N{uJx{q)) + I) cos{n'^^,^T) + N{LJx{q)) cos{n'^ ,^T))], 

phonons, f\{q) ^ {c\q)- The phonon spectrum for ID (34) 
and 2D crystals is plotted in Fig. |8] and Fig.|9] 



Appendix B 

1. Decay of a single excitation 

We start by deriving the decoherence rate 7ip/i for sin- 
gle excitations. To first order in perturbation theory we 



where N{ujx{q)) = \ j (^ehu^(q)/{kBT) _ jg ^j^^ thermal 
occupation number for phonons with the phonon spec- 
trum uj\{q) (see Fig.jsjand Fig.[9| and fi^ ^ — LjJx{q)^i^k- 
LJk is the dispersion relation given in Eq. ( 14 1 . Here we 
have used the property that 
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FIG. 8: Phonon excitation spectrum in ID. a is the separa- 
tion between the molecules. The units of u){q) are Udd/^fP- 





FIG. 9: Phonon Excitation spectrum in 2D for the two dif- 
ferent branches of the phonon spectrum, qx and qy label the 
2D quasimomenta. The units of uo\.2[qx,qy) are Udd/\/13. 



1 ^ 



q \=1 



3a (aA,qe^(''"''-"^(''»* + a^e-'(''-'-----('^»*) 



(35) 



where in the 2D case the vectors e\ are the two orthonor- Ck qX 
mal polarization vectors of the two phonon branches. 



phon~po \'4^k 



and can be explicitly writ- 
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ten as 



- -3, 



2NVf3fxiq) 



E 



a^(e,.(r°-r^")) 
4|rP-r9|5 



(36) 



With the substitution of 
9 



sin(q • r'^) 



r (37) 



Eq. ( 36 I can be simplified to: 



(38) 



In the regime where the Fermi Golden Rule is expected to 
be valid the decay rate is constant; 7(t) = 7 and following 
a similar procedure described in ref. |16j one can show it 
is given by: 



liph 



7r(C + 4i3o)^ 
d^{ak) 



E 



(27r) 



(39) 



gx{k)[{N{LOxik)) + l)<S(wA(fc) - iOk) 



The resonance condition is defined as {LUx{qQ)±uj^±) = 

and assuming that qi]^ = is the only possible solution, 
the decay rate is determined by the fc ^ limit of the 



integrant in Eq. (39 I which is given in ID by: 



7^^ 



(40) 



Due to the finite value of 7^^, and the linear dependence 
on N of ^Trthe probability of remaining on the symmetric 
manifold decreases exponentially with N and decoher- 
ence due to phonons is certainly a limiting factor. 
On the other hand in 2D, the decay rate vanishes as 



7fp"^(X lim(^ + 4So)^A:BTg^0 



(41) 



This conclusion however is only a rough estimation and 
our numerical simulations shows that phonons can in- 
duce important decoherence effects in 2D even in finite 
crystals. 



2. Decay of two dipolar excitations 

In the main body of our paper we have stated that 
72p/i ~ '^"fiph- A detailed derivation is given below. 
To first order in perturbation theory we can write 



-I2ph(t)/h I _ 



I dt' I dr |5fc,fc^,aP[(iV(c.A(g)) + l)cos(^]+^^^,r) + Ar(c^A(q))cos(^]-fc ,,,r 

-^0 q,A,fe,fc'>0 



where ^ ^.i = ^^a('7) i {uik + ujk') and Sk.k'.q.x is rewritten in the same form as 



Sk.k'.q.X — ^34 



3ii 



1 - a4(eA • (rO - r?)) 



2N^fx{q) 



4|rP -r9|5 



2N ^fxiq) ^ 4go)('^q,-fc4',0 + Sq,^k'SkM) - ^^-q.fc + fc') X 



In the above expression the term proportional to ^/iV is much smaller than the one proportional to ^ + ABq and 
can be neglected. Under this assumption, j2ph{t) l2ph-, 



l2ph ~ 2jiph 



(44) 
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